Múltiple return times in the quantum baker map 
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For chaotic classical systems, the distribution of return times to a small region of phase space is universal. We 
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1. Introduction 

In classical dynamical systems, one may study 
how long it takes a system to return to the phase 
space region it occupied at time zero. In gen- 
eral, this time will increase when the phase space 
area under consideration is decreased. M. Kac 
has shown that, for an ergodic dynamical system, 
the average return time, assuming that the sys- 
tem starts in a set A, is given by (x{A)~ l where 
the average is taken over the ergodic measure \x. 
Related to this are results about asymptotically 
rare events. Consider a sequence Ak of subsets 
of phase space such that lim^oo n(Ak) = 0. If 
the time for a point x in phase space to visit Ak 
is denoted by Tk(x), then the question arises how 
the product /j,(Ak)Tk(x) behaves. This turns out, 
for a large class of models, to have an exponential 
distribution ||. 

One may wonder whether a similar situation oc- 
curs in quantum mechanics. In finite dimensional 
systems, looking with increasing precision (corre- 
sponding to the k — > oo límit of above) means 
increasing the dimension of the system, which 
amounts to taking its classical limit. Here, the 
inverse dimension plays the ròle of Planck's con- 
stant. Return times were studied for various 
quantum dynamical systems. A natural way of 
measuring returns is provided by the notion of 
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transition probability: for any two vectors tp and 
tp in a Hilbert space it is defined as |(y|'0)| 2 · 

2. Múltiple return times 

If a system starts from an initial state ífo , one can 
consider 
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where ip n is the state of the system after n time 
steps. For finite dimensional systems this func- 
tion is almost periòdic. Interesting behaviour is 
displayed only when looking at the correct time 
scale. It is known that, for classically chaotic 
systems, the quantum system follows its classical 
counterpart in the sense of coherent states up to a 
certain breaking time of the order of log N, where 
N is the dimension of its Hilbert space of states. 
We propose to analyse the function ([!]) through 
the spectrum of the Gram matrix associated with 
a sequence of normalised vectors ipo, ■ ■ ■ , tpK-i- 
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Gk is a positive semi-defimte matrix and its 
eigenvalues are independent of the order of the 
vectors and of any extra phase added to any of 
thcm. 

Information about the amount and the frequency 
of vectors that return to a small neighbourhood 
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in Hilbert space lies encoded in the spectrum of 
this matrix. This fact is best appreciated by 
examining a simple example. Consider the se- 
quence {ipi,ip2,ip2,4>i,4>3,4>4,4>i}, where the ipi's 
are normalised and mutually orthogonal. This 
amounts to specifying a sequence of letters be- 
longing to a given alphabet (the classical case). 
The Gram matrix consists of only zeroes and ones 
(because of the orthonormality) and its spectrum 
is {0,0, 0,1, 1,2, 3}. An eigenvalue diffcrcnt from 
zero indicates the number of times a certain vec- 
tor is repeated and its multiplicity telis how many 
vectors are repeated that often. In our example, 
only two vectors are repeated, the vector ip\ twice 
and tp2 once, summing up to thrce, which is ex- 
actly the multiplicity of zero. This example pro- 
vides the link with the return times mentioned 
in the introduction: one considers a partition of 
phase space and keeps track of the different sets 
a particle visits during its evolution. The asso- 
ciated Gram matrix has a one for its ijth entry 
if the visited sets at times i and j are the same 
and a zero otherwise. Hence calculating the spec- 
trum of the Gram matrix immediately gives all 
the return times. 

In general, for a true quantum system, the vectors 
in a sequence will not be orthogonal and the over- 
laps will be complex numbers with absolute value 
smaller than or equal to one. However, a similar 
interpretation of the spectrum can be obtaincd 
remembering that overlaps with absolute value 
close to 1 indicate vectors close to each other, 
whilc overlaps close to zero indicate almost or- 
thogonal vectors. Generically, the non-zero eigen- 
values of the Gram matrix are non-degenerate 
and the notion of multiplicity of above should be 
replaced by the spectral density, i.e. a degenerate 
eigenvalue is now replaced by different eigenval- 
ues close to each other. The multiplicity of zero is 
equal to the number of linear dependències in the 
given sequence. Two rather extreme cases will 
clarify this. A large concentration of eigenvalues 
around 1 indicates that a large amount of vectors 
have small mutual overlaps. This can be com- 
pared with the classical situation of above, with 
a spectrum consisting mainly of ones and only 
a few zeroes, indicating very few repetitions. A 
presence of many eigenvalues close to zero sug- 



gests that many vectors lie close to one another, 
leading to many large overlaps. This compares 
to a classical situation where some letters are re- 
peated many times, giving rise to many zeroes in 
the spectrum and only a few large integers. 
A point that has not been stressed yet is the num- 
ber of time steps to consider. The spectrum of 
the Gram matrix will surely contain an eigen- 
value zero when the number K of vectors exceeds 
the dimension N of the Hilbert space. Including 
too many time steps causes a large degeneracy of 
zero. We shall demonstrate, using numerical sim- 
ulations, that also K has to be taken large such 
that the ratio K/N is a fixed number r: time mea- 
sured on the Heisenberg seale. The same sealing 
is relevant for a random model that can be treated 
analytically. 

A natural choice for the initial state is a coherent 
state centred around a point in classical phase 
space. There are several reasons for this. Tak- 
ing the large N limit immediately implies the 
comparison of quantum systems with different di- 
mensions. This also entails that initial states for 
different dimensions have to be comparable. Co- 
herent states exist for every N, occupy a volume 
1 /N of phase space and have a clear limit, namely 
a single point in phase space. Also, a coherent 
state is typically not a eigenvector of the evo- 
lution which would lead to a completely trivial 
spectrum of the Gram matrix. 
The precise object we shall investigate is the em- 
pirical eigenvalue distribution of the Gram ma- 
trix. If a matrix has eigenvalues Ai, • • • , \k, then 
its empirical eigenvalue distribution is 

K 

p(dA) = ^ <5(A - Ai) dA. 

i=l 

In this way, the number of eigenvalues in a set 
Aclis given by p(A). 

We shall illustrate this approach by example. We 
consider two simple quantum dynamical systems 
frcqucntly encountered in the quantum chaos lit- 
erature. They have clear ly defined classical lím- 
its. The first example is the quantum baker map 
and the second the kicked top. 
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3. Models 

The quantum baker map is a dynamics on C N 
with N even m. Defining the discrete Fourier 
transform on í? 1 
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the unitary evolution operator of the baker map 
is 

°-»(? 4)- 

Its classical limit can be shown to be 

</\ / 2q mod 1 \ 
p) " {p/2+[2q]/2j 

where is a point in the phase space [0, 1] x 

[0, 1]. The position and momentum q and p are 
the classical limits of 
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The baker map is an example of a chaotic sys- 
tem: it has a positive Lyapunov exponent equal 
to log2. 

First, we numerically construct a coherent 
state in an iV-dimensional Hilbert space centred 
around a chosen classical point in phase space. 
One way of doing this is by calculating the ground 
state of an appropriate Harper operator and sub- 
sequently shifting it over to the desired classical 
point. Then, we apply the operator U a num- 
ber of times to calculate the overlaps needed in 
the Gram matrix. Note that when the first row 
of the Gram matrix is calculated, all other ma- 
trix elements are known because of the unitarity 
of the evolution. Unitarity has another interest- 
ing consequence: the precision of the calculation 
does not exponentially decrease with the num- 
ber of time steps. One way to see this, is apply- 
ing the unitary evolution to an initial state many 
times, followed by apply ing an equal number of 
times the adjoint of the evolution. The resulting 



state is remarkably close to the original state, de- 
spite the presence of numerical round-off errors 
in the components of the vectors and the unitary 
matrix. This is in sharp contrast with numeri- 
cal simulations of chaotic classical systems where 
positive Lyapunov exponents cause the numeri- 
cal errors to quickly exceed the desired precision. 
The eigenvalues of the Gram matrix are calcu- 
lated and the empirical eigenvalue distribution is 
represented by the bar diagrams of Figure @. The 
sòlid black line will be commented on later in this 
letter. 

The first thing that Figure [ï] reveals is the time 
scale N on which a non-trivial eigenvalue density 
is present. For a fixed time on the Heisenberg 
time scale, i.e. a fixed ratio r = K/N , the eigen- 
value densities are very similar. For comparison, 
choosing K of the order of log N produces a Gram 
matrix close to the identity matrix and therefore 
a spectrum which is greatly concentrated around 
one. For times much larger than N the linear de- 
pendències dominate the picture and the empir- 
ical eigenvalue distribution tends to a Dirac dis- 
tribution at zero. Secondly, we note that we did 
not specify which coherent state we began with. 
Our simulations show that the initial state does 
not affect the form of the eigenvalue distribution. 
Diagrams, very similar to those of Figure were 
obtained for different initial states. 
The same procedure was carried out for a sec- 
ond example: the quantum kicked top. This is 
a quantum spin-j system with bàsic observables 
J x , J y and J z obeying the usual commutation re- 
lations [J x , J y ] = iJ z + cycl. perm.. The unitary 
evolution is given by 



U 



e -ikJ z /2j e -ipJ y 



The classical limit of this system is a dynamics 
on the two dimensional sphere given by 

x = (x cosp + z sinp) cos(fcz cosp — kx sinp) 

+ y sin(fcz cosp — kx sinp), 
y = (xcosp + z smp) sin(fcz cosp — kx sinp) 

— y cos(fcz cosp — kx sinp), 
z = —x sinp + z cosp, 

with (x, y, z) a unit vector in R 3 which is the clas- 
sical limit of (J x , J y , J z )/j- This map is a com- 
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Figure 1. Eigenvalue density of Gram matrices for the quantum baker map for various vàlues of N and 
t = K/N 



position of a rotation around the y-axis over an 
angle p and a torsion around the z-axis over an 
angle proportional to the z-component. Numer- 
ical studies suggest the following behaviour. For 
small vàlues of the parameters k and p, this sys- 
tem displays regular behaviour in large parts of 
the sphere, separated by chaotic regions. The reg- 
ular islands shrink as the parameters grow, lead- 
ing eventually to almost global chaos (See M). 
This is an interesting feature of this model: vary- 
ing the control parameters k and p changes the 
behaviour of the classical model considerably. It 
can be concluded from Figure || that this differ- 
ence in behaviour is also present in the quantum 
kicked top. For small vàlues of the control param- 
eters k and p there are a lot of approximate re- 
turns during the first K — tN time steps. This is 
clear from the large number of eigenvalues close to 
zero. We expect that in the limit K = tN — > oo 
this distribution tends to a degenerate one in zero 
as the system visits only a small fraction of the 
available states. Enlarging the control parame- 



ters places more weight on higher eigenvalues and 
a totally different eigenvalue distribution appears. 
These distributions appear to be very similar to 
the ones for the baker map for the same value of 

T. 

One unexplained feature in the diagrams is the 
sòlid line. This is actually the graph of an ana- 
lytical function obtained in the following way. In- 
stead of trying to calculate the limiting empirical 
eigenvalue distribution for a true quantum me- 
chanical model, like the baker map or the kicked 
top, one can try to simplify the calculations by 
considering a random vector model. For an N- 
dimensional Hilbert space, a sequence of K nor- 
malised vectors is chosen randomly and indepen- 
dently according to the uniform distribution on 
the unit sphere. In this way, one can try to mimic 
the evolution of a system which wanders wildly in 
Hilbert space. The limiting eigenvalue distribu- 
tion (N,K — > oo and K/N -> r) of the Gram 
matrix associated with this sequence of random 
vectors can be calculated explicitly and is given 
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Figure 2. Eigenvalue density of Gram matrices for the quantum kicked top for various vàlues of r, k en 
p (r is defined as the ratio K/N with N = 2 j + 1 = 201) 



by the Marchenko-Pastur distribution 
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This density is the sòlid line depicted in Figures ^ 
and |[ A possible explanation for the similar- 
ity between the eigenvalue distributions gener- 
ated by the baker map and the kicked top and 
the Marchenko-Pastur distribution is that, whcn 
the number of time steps is linear in the dimen- 
sion of the Hilbert space, subsequent vectors in 
a chaotic quantum evolution are almost random 
and visit all of the Hilbert space without prefer- 
ence for one direction or region. 

4. Conclusion 

The eigenvalue distribution of the Gram matrix 
appears to be a meaningful object in the so-called 
deep quantum regime where the system is allowed 
to visit a finite fraction of the total number of 



available states. Moreover, for systems with a 
chaotic classical limit, numerical simulations in- 
dicate that the limiting distribution is universal 
and coincides with that of random vectors, even 
if the latter model is not a unitary dynamics. If 
the classical limit is regular, we find a degener- 
ate distribution. Linking the limiting spectrum of 
the Gram matrix with the distribution of energy 
or quasi-energy level separations definitely de- 
serves further research. Another interesting ques- 
tion, that is also relevant in quantum information 
theory, is what happens when the dynamical sys- 
tem is no longer isolated but weakly coupled to 
an environment. In such a situation, density ma- 
trices instead of vector states must be used to 
describe the decoherence induced by the environ- 
ment. 
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